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ABSTRACT 

We present stellar kinematics and orbit superposition models for the central regions of four Brightest 
Cluster Galaxies (BCGs), based upon integral-field spectroscopy at Gemini, Keck, and McDonald 
Observatories. Our integral-field data span radii from < 100 pc to tens of kiloparsecs, comparable to 
the effective radius of each galaxy. We report black hole masses, M,, of 2.lli['g x 10^° Mq for NGC 
4889, Q.lt^l X IO^Mq for NGC 3842, and 1.3l[J;^ x 10^ Mg for NGC 7768, with errors representing 
68% confidence limits. For NGC 2832 we report an upper limit of M, < 9.0 x 10^ Mq. Our models of 
each galaxy include a dark matter halo, and we have tested the dependence of M, on the model dark 
matter profile. Stellar orbits near the center of each galaxy are tangentially biased, on comparable 
spatial scales to the galaxies' photometric cores. We find possible photometric and kinematic 
evidence for an eccentric torus of stars in NGC 4889, with a radius of nearly 1 kpc. We compare our 
measurements of M, to the predicted black hole masses from various fits to the relations between 
M, and stellar velocity dispersion (u), luminosity (L), or stellar mass (M^,). Still, the black holes 
in NGC 4889 and NGC 3842 are significantly more massive than all cr-based predictions and most 
L-based predictions. The black hole in NGC 7768 is consistent with a broader range of predictions. 



ije.g., 'Bolton' 1972': 'Schodcl et al.' '2002'; 
3^ [Bender ct al. 2005; Hcrrnstcin ct al. 



1. INTRODUCTION 

For four decades, dynamical studies have provided ev- 
idence for the existence of black holes. In the most 
exquisite cases, the orbital motions of individual stars 
or megamasers place strong lower limits on the density 
of a central object an d rule out virtually all alternatives 
to a black hole 
Chez et all [200! 
2005[) . More generally, massive dark objects have been 
dynamically detected in a rapidly growing number of 
galactic nuclei and are widely assumed to be black holes 
(e.g., Sargent et al. 1978; Tonry 1987; for reviews and 
compilations see Kormendy & Richstone 1995; Giiltekin 
et al 2009a; McConnell et al. 2011b). The menagerie 
of black holes spans from stellar-mass objects to "super- 
massive" behemoths, whose masses, M,, can approach 
10^° Mq. Strong dynamical evidence for "intermediate- 
mass" black holes with M, ^ 10'^ — 10^ Mq is more elu- 
sive but has been recorded in a few stellar systems (Geb- 
hardt et al. 2005; Liitzgendorf et al. 2011; Jalali et 
al. 2012; c.f. Baumgardt et al. 2003; van der Marel & 
Anderson 2010). 

In galactic nuclei, Af, appears to correlate with 
the stellar mass, luminosity, and velocity disper- 
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tier goal is to explore the extrema of these black 
hole scaling relations by directly measuring M, in 
low-mass galaxies and in extremely massive galaxies. 
Here we focus on Brightest Cluster Galaxies (BCGs), 
which are among the most massive galaxies in the 
present-day Universe. BCGs typically reside deep in 
the gravitational potentials of rich clusters, and their 
environment may provide a unique path for galaxy 
and black hole growth. BCGs have been observed 
to follow a steeper relation between luminosity (L) 
and velocity disp ersion {a) than less massive ellipti- 
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Bernar di et al. '2007; 'Desroches et alj 
2007a : von der Linden ct al. 2007). 



Bernardi etal.1 (|2007f )~and iLauer et all (|2007al ) have 
noted that the M, — L relation predicts systematically 
more massive black holes than M, — a for the most mas- 
sive galaxies. 

To date, there are eight groups or clusters where M, 
has been measured dynamically in the massive central 
galaxy: Coma (NGC 4889; McConnell et al. 2011b), 
Fornax (NGC 1399; e.g. Houghton et al. 2006; Geb- 
hardt et al. 2007), Virgo (M87; Gebhardt et al. 2011), 
Abell 1367 (NGC 3842; McConnell et al. 2011b), Abell 
1836 (PGC 49940; Dalla Bonta et al. 2009), Abell 2162 
(NGC 6086; McConnell et al. 2011a), Abell 3565 (IC 
4296; Dalla Bonta et al. 2009) and the IC 1459 group 
(IC 1459; Cappellari et al. 2002). Only three of these 
systems (Coma, Abell 1367, Abell 3565) are rich clusters. 
In Fornax and Virgo, the central cD galaxy is not even 
the brightest member. The brightest galaxy in Fornax 
is NGC 1316, which lies near the cluster outskirts and 



TABLE 1 
Global properties of target galaxies 



Galaxy 


D 

(Mpc) 


Mv 


rinf 

(") 


(") 


0"off 

(kms-1) 


(°) 


b/a 


{°) 


Features 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 



NGC 4889 


103.2 


-23.91 


1.5 


54.5 


347 


80 


0.73 


90 




NGC 3842 


98.4 


-23.21 


1.2 


37.8 


270 


5 


0.86 


90 


Isophotal twist at r > 12" 


NGC 7768 


112.8 


-22.92 


0.14 


23.1 


257 


60 


0.75 


90 


Nuclear dust disk 


NGC 2832 


101.9 


-23.87 


< 0.7 


67.4 


334 


-21 


0.68 


90 


Isophotal twist at r < 8" 



Notes: Column 1: galaxy. Column 2: co-moving distance, derived from cluster systemic velocity with Ho — 70kms ^ Mpc ^. 
Column 3: absolute l/-band magnitude. Column 4; black hole radius of inf luence, derived from our best-fit value of Mm and 
the effective velocity dispersion in column 6. Column 5: effective radius from lLauer et all (|2007al ). Column 6: effective velocity 
dispersion, from luminosity- weighted measurements between rinf and r^g. Column 7: major-axis position angle, measured from 
north toward east. This angle is used to define the equatorial plane of each galaxy in our stellar orbit models. Column 8: 
observed axis ratio. Column 9: inclination assumed for stellar orbit models. 90° corresponds to edge-on inclination. Column 
10: notable photometric features. 



hosts a black hole with M, = 1.5 x 10^ Mq (jNowak et alj 
|2008[ ). In Virgo, M49 anchors a sub-group more than 1 
Mpc from the more centralized M87 and M86 sub-groups, 
and hosts a black hole with M, = 1.5 x 10^ Mq (Shen 
et ah, in prep.). In order to thoroughly explore black 
hole and galaxy co-evolution in different environments, 
we must measure M, in a larger sample of BCGs, with 
more examples from rich galaxy clusters. This paper 
includes expanded discussion of the black hole measure- 
ments in the BCGs of Coma and Abell 1367. 

This paper marks our continued effort to measure M, 
using stellar dynamics in BCGs beyond the Virgo cluster. 
BCGs are rare objects and typically lie at large distances, 
making high-resolution observations difhcult. Addition- 
ally, their centers are typically fainter than less massive 
ellipticals, and 8- to 10-meter telescopes are required to 
obtain high-quality spectra at angular scales comparable 
to the black hole radius of influence, rjnf = GM,a^^. 
For targets at distances ^ 100 Mpc and predicted black 
hole masses ^ 10^ Mq from the M, — a relation, adaptive 
optics (AO) is necessary to resolve rinf '^ 0.1". Under 
good conditions, seeing-limited observations can resolve 
the gravitational influence of extremely massive black 
holes (M, ~ IQ^^ AIq). Wide-field kinematic measure- 
ments are necessary to trace the galaxies' stellar mass 
profiles and dark matter halos. 

We report measurements of M, and the i?-band stellar 
mass-to-light ratio, M*/L/j, in four BCGs at distances 
of - 100 Mpc: NGC 4889 of the Coma cluster (Abell 
1656), NGC 3842 (AbeU 1367), NGC 7768 (Abell 2666), 
and NGC 2832 (AbeU 779). We have obtained high- 
resolution data of the line-of-sight stellar velocities in 
the central regions of the four galaxies using instruments 
on the Gemini North and Keck Telescopes. In addi- 
tion, we have used the 2.7-meter telescope at McDon- 
ald Observatory to measure stellar kinematics at large 
radii. At all spatial scales, we use integral-field spectro- 
graphs (IFSs) to obtain full 2-dimensional spatial cov- 
erage, which places tighter constraints on stellar orbits. 
We determine M, and Mi,/ Lb, with axisymmetric orbit 
superposition models, which include a dark matter com- 
ponent in the gravitational potential. We have reported 
the bla ck hole measurements o f NGC 4889 and NGC 
3842 in lMcConnehet al.l (|2011bD . but due to space hmi- 



tations, only the basic information was presented there. 
Here we provide a comprehensive discussion of the data 
analysis procedures, kinematic information, and stellar 
orbit modeling for all four galaxies. 

The paper is organized as follows. In §2, we describe 
our photometric and spectroscopic observations and data 
reduction procedures. In §3, we describe our proce- 
dures for extracting two-dimensional kinematics from 
the data obtained on three IFSs (CMOS, OSIRIS, and 
the Mitchell Spectrograph). The outputs of this step 
are non-parametric line-of-sight velocity distributions 
(LOSVDs) in two-dimensional spatial bins in the central 
regions of each galaxy. In §4, we show maps of the lowest 
four Gauss-Hermite moments of the LOSVDs for the four 
galaxies. We discuss both the two-dimensional maps and 
radial profiles of the kinematic moments. In §5 we sum- 
marize the stellar orbit modeling procedure. We report 
our measurements of Af, and Mi,/ Lji and describe how 
these measurements depend on the assumed dark matter 
halo profile, as well as other possible systematic biases. 
We also discuss the relative contributions of radial and 
tangential orbits in each galaxy. In §6 we summarize our 
results and compare our measurements of M, to predic- 
tions from the Af, — a and Af, — L relations. Throughout 
this paper, we assume i?o = 70 km s^^, Vlrn ~ 0.27, and 
r^A = 0.73. 

2. OBSERVATIONS 
2.1. Global Galaxy Properties 

We list the basic properties of each galaxy in Table [TJ 
We compute co-moving distances from the average veloc- 
ity of the host galaxy cluster with respect to the cosmic 
microwave backgro und (CMB). We us e helio centric clus- 
ter velocities from iLauer fc PostmanI ()1994| ) and trans- 
late them to the CMB rest frame with the NASA/IPAC 
Extragalactic Database (NED). To compute the V- 
band luminosity, Ly, for NGC 4889, NGC 3842, and 
NGC 2832, we adj ust the absolute magnitudes from 
ILauer et al.l (j2007aU bf) to our assumed distances. For 
NGC 7768, we adopt the ap parent y-band magnitude o f 
12.40 from the RC3 catalog (|de Vaucouleurs et al.|[T99l 
and assume a luminosity distance of 115.8 Mpc. 

Galaxies with spatially resolved stellar kinematics are 
often placed on the Af, — a relation using OeS, the 
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Fig. 1. — De-projected ij-band stellar luminosity density versus radius along the major axis (solid line) and minor axis (dotted line) of 
each galaxy. The dashed vertical lines mark the outermost extents of photometric data from HST and KPNO. Luminosity densities beyond 
115" are derived from a de Vaucouleurs surface brightness profile, (a) NGC 4889. (b) NGC 3842. (c) NGC 7768. (d) NGC 2832. 



luminosity-weighted average velocity dispersion out to 
one effective radius (roff). Yet when the black hole ra- 
dius of influence, rinf , subtends a large angle, the con- 
ventional measurement of Ceff can include a direct signal 
from the black hole. In extreme cases, the velocity dis- 
persion only reaches its "average" value at r < rinf (see, 
e.g., Gebhardt et al. 2011 for M87). As the M, - a 
relation is widely interpreted as an empirical correlation 
between independent parameters, it is important to re- 
move the black hole's direct influence on (Toff. For our 
BCGs, we measure the luminosity-weighted average ve- 
locity dispersion for Hnf <r< reff : 



2 



/ro 



'^^ + «rld) Hr)dr 



/:: I{r)dr 



(1) 



where /(r) is the galaxy's one-dimensional stellar surface 
brightness proflle. 

We compute r-mf and Ucfi from the integral-fleld kine- 
matics and black hole masses presented herein; at each 
radius, we average Wj-ad and a over all polar angles sam- 
pled. Because rinf = GM,a~^, we compute ri„f and CTcff 
iteratively; the iterations converge quickly in all cases. 
For NGC 2832, we only have upper limits for M, and 
Hnf, and we assume nnf = 0. For NGC 4889 and NGC 
3842, our upper integration limit is the maximum radius 
of our kinematic data, which corresponds to 0.4 roff and 
0.8 Toff, respectively. 

2.2. Photometry 

We measure the stellar light profiles of each BCG with 
a combination of i?-band (0.6 fim) and /-band (0.8 fim) 
photometry. For radii, r, out to 10" we adopt high- 



resolution surface brightness profile s from WFP C2 on 
the Hubble Space Telescope (HST) ()Laine et all [2003). 
At larger radii we use _R-band data from Lauer, Post- 
man & Strauss (private communication), obtained with 
the 2.1-m telescope at Kitt Peak National Observatory 
(KPNO). The KPNO data have a field-of-view (FOV) of 
5.2' X 5.2', which enables accurate sky subtraction. Our 
ground-based surface brightness profile of each galaxy 
extends to 115". We have combined the individual pro- 
files from WFPC2 and KPNO at overlapping radii be- 
tween 5" and 10", varying the respective weights such 
that the WFPC2 data contribute 100% to the combined 
profile at r = 5" and the KPNO data contribute 100% 
at r = 10". Before stitching the profiles, we translate 
the WFPC2 profile to i?-band using the average R — I 
color bet ween 5" and 10" NGC 7768 has a dust disk at 
r < 0.5" (jGrillmair et al.|[l994D . so we measure the in- 
nermost stellar profile with additional iJ-band (1.6 piva) 
photometry from /T^T/NICMOS (1997; PI Tonry). 

We convert surface brightness to stellar luminosity den- 
sity u sing the deprojection procedure of IGebhardt et al.] 
(jl996[) . In order to model each galaxy, we must define a 
symmetry axis. The symmetry axis in three dimensions 
is projected to the minor axis on the sky, and the model 
galaxy's equatorial plane corresponds to the photomet- 
ric major axis. Isophotal twists are evident in the outer 
part of NGC 3842 and the inner part of NGC 2832; in 
each case, we choose the major axis position angle that 
matches the largest range of radii covered by our kine- 
matic data. The major axis of NGC 4889 is near 80° 
east of north at all radii. We summarize our photomet- 
ric measurements of each galaxy in Table [1] The major- 
and minor-axis luminosity density profiles are shown in 



TABLE 2 
Spectroscopic observations 



Galaxy 


Instrument 


Date 


tint 


A Range 


PA 


Ax 


FWHM 








(S) 


(nm) 


{°) 


{") 


(") 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




NGC 4889 


OSIRIS 


May 2008 


9x 900 


1473 - 1803 


135 


0.05 






NGC 4889 


OSIRIS 


May 2009 


9 X 900 


1473 - 1803 


80 


0.05 






NGC 4889 


CMOS 


March 2003 


5 X 1200 


744 - 948 


90 


0.2 


0.4- 


0.7 


NGC 3842 


OSIRIS 


May 2010 


9 X 900 


1473 - 1803 


10 


0.05 


0.06, 


0.19 


NGC 3842 


CMOS 


April 2003 


5 X 1200 


744 - 948 


-30 


0.2 


0.4- 


0.7 


NGC 3842 


Mitchell 


March 2011 


3 X 1200 


358 - 588 





4.1 






NGC 7768 


OSIRIS 


September 2010 


14 X 900 


1473-1803 


-125 


0.05 


0.04, 


0.09 


NGC 7768 


Mitchell 


September 2011 


13 X 1200 


358 - 588 





4.1 






NGC 2832 


CMOS 


March 2003 


4 X 1200 


744 - 948 


30 


0.2 


0.4- 


0.7 


NGC 2832 


Mitchell 


January 2008 


11 X 1200 


354 - 584 





4.1 






NGC 2832 


Mitchell 


February 2008 


28 X 1200 


354 - 584 





4.1 







Notes: Column 1: galaxy. Column 2: instrument. OSIRIS (OH-Suppressing Infra-Red Imaging Spectrograph) was used on 
Keck II with LGS-AO. CMOS (Gemini Multi-Object Spectrograph) was used on Gemini North. The Mitchell Spectrograph 
was used on the McDonald Observatory 2.7-m telescope. Column 3: observing dates. Column 4: number of science exposures 
X integration time per exposure. Column 5: observed wavelength range. Column 6: position angle of the long axis for OSIRIS 
and GMOS, measured from north toward east. Column 7: angular size of lenslets or fibers. Column 8: PSF FWHM at science 
wavelengths. For OSIRIS observations of NGC 3842 and NGC 7768, we list the FWHM values of the inner and outer PSF 
components. 



Figure [TJ 

We use i?-band photometry in our stellar orbit mod- 
els and therefore constrain Mi,lLii. It is useful to de- 
rive Mi^jLy so we can compare our results with other 
studies of mass-to-light ratios in early- type galaxies. 
To c onvert from R- to F-ba nd, we use galaxy colors 
from lPostman fc Laued (|1995[ ) and the Sloan Digital Sky 
Survey , with filter translations from Blanto n & Rowcis 
(|20(T7fl . 

2.3. Integral-field Spectroscopy 

We measured kinematics in NGC 3842, NGC 4889, 
NGC 7768, and NGC 2832 using three different 
integral-field spectrographs (IFSs) to cover two or- 
ders o f magnitude in radius. Data from GM OS- 
North (|Allington-Smith et al.ll2002t iHook et al.ll2004[ ) on 
the 8-m Gemini Observatory North telescope provide 
the most reliable measurements of the central stellar 
kinematics in each BC G. Additional observations with 
OSIRIS (Larkin et al.l 12006 ) and the laser guide star 
adaptive optics (LGS-A O) system (jvan Dam et al.|[2006l : 
iWizinowich et al.ll2006t) on the 10-m W. M. Keck II tele- 
scope provide higher spatial resolution but substantially 
worse signal-to-noise. We measured wide-field kinemat- 
ics with the George and Cyn thia Mitchell Spectrograph 
(formerly VIRUS-P; iHill et al. 2008) on the 2.7-m Har- 
lan J. Smith Telescope at McDonald Observatory. Our 
spectroscopic observations are summarized in Table [H 

2.3.1. GMOS 

GMOS-N is a multi-purpose spectrograph on Gemini 
North. We used GMOS in IPS mode, which maps the 
science field and a simultaneous sky field with hexagonal 
lenslets. We observed the center of each BCG with the 
CaT filter, centered near the infrared Call triplet. A 



representative spectrum from each galaxy is shown in 
Pigure [21 The GMOS data were reduced using version 
1.4 of the Gemini IRAP software packag^. 

GMOS is seeing-limitcd, with a lenslet scale of 0.2". 
We measure full widths at half-maximum (PWHM) of 
0.4 — 0.7" for point sources in acquisition frames from 
the GMOS imager. Por all stellar orbit models discussed 
below, we have assumed seeing of 0.4" FWHM for GMOS 
data. In a preliminary trial with NGC 3842, we obtained 
consistent values of M, and Af*/L from stellar orbit mod- 
els assuming 0.4" and 0.7" seeing. 

2.3.2. OSIRIS 

OSIRIS is a near- infrared (NIR), IPS built for use with 
the Keck AO system. We observed each BCG with the 
0.05" lenslet scale and broad 7J-band filter, which cov- 
ers several metal absorption lines and CO and OH vi- 
brational bandheads at observed wavelengths of 1.47- 
1.80 /im. Pigure [3] shows representative OSIRIS spec- 
tra. OSIRIS has no sky lenslets, and extended objects 
such as BCGs cover the entire science field. We observed 
each BCG with a repeated "object-sky-object," dither 
sequence, such that every 900-s science frame was imme- 
diately preceded or followed by a 900-s sky frame. We 
observed kinematic template stars using the same filter 
and lenslet scale as the BCGs. Spatial variations in in- 
strumental resolution are negligible relative to the veloc- 
ity broadening in BCGs. 

We used version 2.3 of the OSIRIS data reduction 
pipelineQto perform sky subtraction, spatial flat-fielding, 
spectral extraction, wavelength calibration, and spatial 

^ available from Gemini Observatory, at 

|http://w ww.gemini.edu/sciops/data-and-results/proccssing-s oftware| 

' aTOilable from the UC LA Infrared Laboratory, at 
|http://ir lab. astro. ucla.edu/osiris/pipeline. html 



mosaicking of 3-dimensional data cubes. We used custom 
routines to remove bad pixels and cosmic rays from 900- 
s exposures and to calibrate for telluric absorption. We 
estimated the AO PSF by observing the LGS-AO tip/tilt 
star with the OSIRIS Imager. For NGC 3842 and NGC 
4889, we typically measured the PSF once per half-night. 
For NGC 7768 we measured the PSF four times over an 
observing span of eight hours. The weighted average PSF 
for NGC 7768 has a Strehl ratio of 22% and is well-fit 
by an inner component with 0.04" FWHM and an outer 
component with 0.09" FWHM. In a similar i nvestig ation 
of the BCG NGC 6086, iMcConneh eTaH (|2011ari var- 
ied the AO PSF and obtained a systematic error of only 
~ 10% in the fi nal measuremen t of M «. In their investi- 
gation of M87, iGebhardt et al.l (|2011[ ) measured similar 
values of M, for PSFs with different Strehl ratios. 

OSIRIS data for NGC 4889 and NGC 3842 have lower 
quality than we initially expected. Even spectra with 
relatively high signal-to-noise (S/N) yield unconvincing 
kinematic measurements. One severe noise source is tel- 
luric OH emission, which varies faster than the 900- 
s exposure time and leaves large residuals even after 
subtracting sky frames. Although we mask the narrow 
residual features when measuring stellar kinematics, the 
missing spectral channels conceal information about the 
overall shapes of absorption features. Further difhculties 
arise from defining a continuum level among the blend 
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Fig. 2. — Example spectra of the calcium triplet region from the 
GMOS IFS. For each galaxy, the central spectrum (r < 0.25") is 
displayed, as well as spectrum near the major axis, corresponding 
to an inner radius of 1.1" and an outer radius of 1.7". The major- 
axis spectra correspond to the east side of NGC 4889, the north 
side of NGC 3842, and the north side of NGC 2832. The thick 
cyan line for each spectrum represents the best-fitting, LOSVD- 
convolved stellar template. Dotted lines indicate spectral channels 
that were masked during fitting. Each spectrum was continuum- 
divided before fitting and plotting. 



of metallic and molecular absorption features in iJ-band. 
Finally, OSIRIS suffered from elevated detector temper- 
atures throughout 2009, including the bulk of our ob- 
servations of NGC 4889. Data from this period exhibit 
increased dark noise and possible errors in spectral ex- 
traction and wavelength calibration. 

Our OSIRIS data for NGC 4889 do not have suffi- 
cient quality to include in stellar orbit models. We have 
run stellar orbit models of NGC 3842 with and without 
OSIRIS data; we compare results from these models in 
Section [5.2.21 OSIRIS data were included in all models 
of NGC 7768. 

2.3.3. Mitchell Spectrograph 

The Mitchell Spectrograph, formerly named VIRUS- 
P, is an optical IFS on the 2.7-m telescope at McDonald 
Observatory. It features 246 fibers spaced evenly across a 
107" X 107" field of view, with a one-third filling factor. 
Each fiber has a diameter of 4.1"; this is the limiting 
factor in spatial resolution, rather than optical seeing. 

We observed each galaxy with the low-resolution grism, 
which provides wavelength coverage of ~ 360-580 nm, in- 
cluding the Ca H -I- K region, the G-band region, H/3, 
the Mg b region, and several Fe absorption features. 
The spectral resolution varies over different fibers and 
also with wavelength. Our observations using the low- 
resolution grism span instrumental resolution values of 
« 0.45 — 0.65 nm FWHM, with a corresponding range of 
a ^ 100 — 170 km s^^. This is sufficient to resolve the ve- 
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Fig. 3. — Example H-hand spectra from OSIRIS. The central 
spectrum for NGC 7768 covers O.l" X O.l". The central spectrum 
for NGC 3842 corresponds to r < 0.25". We also display an ofl"- 
center spectrum for each galaxy; the corresponding spatial region 
is near the major axis, with an inner radius of 0.4" and an outer 
radius of 0.7". The major-axis spectra correspond to the east side 
of NGC 7768 and the north side of NGC 3842. The thick cyan line 
for each spectrum represents the best-fitting, LOSVD-convolved 
stellar template. Dotted lines indicate spectral channels that were 
masked during fitting, and dashed red lines indicate residual fea- 
tures from telluric OH, which were also masked. Each spectrum 
was continuum-divided before fitting and plotting. OSIRIS spectra 
for NGC 4889 were compromised by warm detector temperatures 
in 2009. 
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Fig. 4. — Example spectra from the Mitchell Spectrograph. For each galaxy, we display spectra corresponding to the smallest and largest 
radii where we employ Mitchell Spectrograph data in our stellar orbit models. The corresponding bins span 3.8 — 5.6" and 24.5 — 35.3" for 
NGC 3842, 3.0 - 4.8" and 18.3 - 28.5" for NGC 7768, and 1.7 - 2.5" and 70 - 101" for NGC 2832. Dotted lines indicate spectral channels 
that were masked during fitting. Each spectrum was continuum-divided before fitting and plotting. 



locity profiles of BCGs {a ^ 250 — 350 km s^^). Because 
we employed the Mitchell Spectrograph for wide-field 
coverage rather than fine spatial resolution, we did not 
perform any sub-dithers to fill the gaps between fibers. 

We used the Vaccine data reduction pipeline to per- 
form bias subtraction and flat-fielding, compute the 
wavelength solution for each fiber, extract a spectrum 
for each fiber, model and subtract the sky spectrum, and 
reject cosmic rays. Detailed descrip tions of the data re - 
duc tion procedures can b e found in lAdams et al.l ()2011[ ) 
and lMurphvet all (|201lD . 

The maximum radius at which Mitchell Spectrograph 
data can yield robust kinematic measurements depends 
on the surface brightness profile of the galaxy, the back- 
ground sky conditions, and the overall integration time 
at a given pointing. For each galaxy, we have binned nu- 
merous fibers at large radii and preserved data for which 
a good kinematic fit was recoverable for at least one bin- 
ning scheme. Following this procedure, our outermost 
bins for kinematic extraction cover radii of 24.5 — 35.3" 
for NGC 3842, 18.3 - 28.5" for NGC 7768, and 70 - 101" 
for NGC 2832. Sample spectra from the central and outer 
regions of each galaxy are shown in Figure |4l 

3. EXTRACTING STELLAR KINEMATICS 

Our dynamical models fit weighted and super- 
posed stellar orbits to line-of-sight velocity distributions 
(LOSVDs) extracted from spectroscopic data. We ex- 
tract LOSVDs with a Maximum Penalized Likelihood 
(MPL) technique, which fits an LOSVD-convolved stel- 
lar template to each galaxy spectrum. The LOSVDs are 



non-parametric, defined at 15 radial velocity bins in our 
orbit mo dels. The MPL fitting method is descr i bed i n 
deta il in iGebhar dt et al.l (|2000br ). iPinknev et all ([200l . 
and iNowak eta l. (2008), and adapta tions for integral 
field da ta of BCGs are described in iMcConnell et al.l 

il20ITa . 

In order to attain sufficient signal-to-noise (S/N) for ef- 
fective kinematic extraction, we perform spatial binning 
on each data set. For CMOS data of NGC 4889, NGC 
3842 and NGC 2832, our central bin combines 7 hexago- 
nal Icnslcts; the corresponding bin diameter is 0.55". For 
OSIRIS data of NGC 7768, we combine 4 square lenslets 
to obtain a central bin of 0.1" x 0.1". At large radii (typ- 
ically r > 15") we combine Mitchell Spectrograph data 
from multiple angular bins. 

A spectral binning factor is necessary to smooth over 
channel-to-channcl noise in some spectra. We typically 
use smoothing factors of 2 spectral pixels for GMOS data 
and 8-10 spectral pixels for OSIRIS data. Lower smooth- 
ing factors yield jagged LOSVDs. Mitchell Spectrograph 
spectra do not require smoothing. 

The near-infrared Call triplet has a clearly defined 
continuum and yields robust kinematic measurements 
across a broad range of stellar templates ( Dressier 1981; 
iBarth et al. 2002 ). Therefore, we fit GMOS spectra with 
a single G9III template star. In spectra of NGC 4889 and 
NGC 2832, the 866 nm calcium line is contaminated by a 
detector artifact and is not included our fits (see Fig. [2]) . 
The 850 nm and 854 nm lines are included in all fits. 

For OSIRIS spectra, we optimize our LOSVD extrac- 
tion by carefully choosing a subset of the many absorp- 



tion features in H-hand. Although we are able to per- 
form a global equivalent width adjustment before fitting 
stellar templates, some of the individual metal lines ex- 
hibit further mismatch between the galaxy and template 
equivalent widths. Some molecular features are blended 
with the suspect metal lines, and others are compromised 
by residuals from telluric OH subtraction. The consis- 
tently reliable features for kinematic extraction are the 
f = 3-6, ly = 4-7, and ly = 5-8 bandheads of ^^CO, span- 
ning the sub-region from 1.61 to 1.68 fim rest, and the 
Mgl absorption feature near 1.504 /im rest (see Fig. [3]). 
We fit OSIRIS spectra using a set of nine template stars 
with spectral types from G8 to M4. Template weights 
are allowed to vary freely and are fit simultaneously with 
the LOSVD. 

Mitchell Spectrograph spectra contain a multitude of 
absorption features, which are dominated by different 
stellar types and require different equiva lent width ad- 
j ustm ents. Following the procedure of iMurphv et al.l 
()2011|) . we divide each spectrum into five spectral sub- 
regions and extract LOSVDs independently for each sub- 
region. Our fits use a set of 16 template stars selected 
from the Indo-US library (Valdes et al. 2004), with spec- 
tral types from B9 to MB. We compute the instrumen- 
tal resolution as a function of wavelength in each fiber 
and convolve the template spectra with an appropriately 
weighted instrumental resolution profile for each galaxy 
spectrum. After fitting each galaxy spectrum, we discard 
sub-regions with visually unconvincing fits or severely 
asymmetric LOSVDs. The LOSVDs from the remaining 
sub-regions are averaged to produce a final representative 
LOSVD. 

We perform Monte-Carlo trials to determine the un- 
certainties in the best-fit LOSVDs. In each trial, ran- 
dom noise is scaled to the root-mcan-squared residual of 
the original fit and added to the galaxy spectrum be- 
fore re-fitting. In each bin of the LOSVD, uncertainties 
are computed from the distribution of trial values. Our 
fitting method sometimes yields LOSVDs with noise at 
large positive or negative velocities. Consequently, we 
force our uncertainties to extend to zero power in the 
LOSVD wings. For MitcheU Spectrograph LOSVDs, we 
adopt the average uncertainty from the individual spec- 
tral sub-regions, or if greater, the l-tr scatter in the sub- 
regions' LOSVD values. 

4. TWO-DIMENSIONAL KINEMATICS 

For each galaxy, we fit a 4th-order Gauss-Hermite poly- 
nomial to each non-parametric LOSVD in order to illus- 
trate the stellar kinematics. Below we highlight spatial 
trends in the Gauss-Hermite moments Wrad, cr, /13, and h^. 
To measure the Gauss-Hermite moments, the LOSVD is 
fit with a function f{v), defined by 
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4.1. NGC 4889 

Figure [S] shows two-dimensional maps of Uradj cTi ^3: 
and /i4, from CMOS observations of NGC 4889. Fig- 
ure[n]also includes radial profiles of each moment in NGC 
4889. The radial p rofiles include our G MOS data and 
measurements from lLoubser et all ()2008[ ). recorded with 
the long-slit spectrograph ISIS at the William Herschel 
Telescope (WHT). 

The kinematic moments in NGC 4889 show several 
asymmetries with respect to the major and minor axis. 
Asymmetries in the line-of-sight velocity dispersion, a, 
likely present the greatest difficulties for accurately mea- 
suring M,. On the east side of the galaxy, a peaks at 
413 ± 22kms^^ near the major axis and remains above 
400 km s^^ through an extended region in the southeast 
quadrant. On the west side, a exceeds 385kms~^ in 
only a single spatial bin, reaching 406 ± 18kms"^ at 
0.6" toward the northwest. At the very center of NGC 
4889, we measure a = 344 ± 16kms~^. A strong deficit 
of radial orbits can produce the central drop in line-of- 
sight velocity dispersion even in the presence of a su- 
permassive black hole; this is further discussed in Sec- 
tion 15.31 Central decreases in velocity dispersion have 
been observed in several other early- type galaxies (e.g. 



van der Marcl'1994':'Pinkncv ct al."2003': 'Houg hton et alj 
2006; Gebhardt et al. 2007; Nowak ct al. 20081). 

The ft.4 moment, which describes whether an LOSVD 
is boxy (/14 < 0) or peaky (ft.4 > 0), varies with radius 
and polar angle in NGC 4889. At radu from 0.7" to 
3.5", hi is negative near the major axis and nearly zero 
toward the minor axis. We have fit and subtracted a 
Gaussian profile from each non-parametric LOSVD and 
have examined the residual velocity profiles. In most of 
the spatial bins near the major axis, the residual profiles 
contain peaks near ±450kms~^, as exemplified in Fig- 
ure [6l An LOSVD sub-component with velocity peaks 
near ±450 km s^^ could correspond to a torus or thick 
ring of stars, including rotating and counter-rotating 
populations. In an eccentric ring, a would exhibit a local 
maximum near periapsis. 

To further investigate the hypothesis that NGC 4889 
hosts a central stellar torus, we have carefully exam- 
ined the inner light profile of the galaxy. We have de- 
termined a best-fit surface brightness profile by fitting 
concentric ellipses to a re duced and de-conv olved image 
from HSTIWYVC2, as in lLaine et all (pOOl . Subtract- 
ing the best-fit ellipses yields a 2-dimensional residual 
image of NGC 4889, displayed in Figure [71 The resid- 
ual image shows two peaks along the major axis, which 
straddle the galaxy center and contain approximately 1% 
of the total light at their corresponding positions. These 
peaks resemble a diffuse torus of stars in the equatorial 
plane of the galaxy, spanning radii from 1" to 2". On the 
east side of NGC 4889, the residual photometric peak 
overlaps with the velocity dispersion peak but does not 
extend to match the extended pattern of high a in the 
southeast quadrant. Boxy LOSVDs overlap with both 
residual photometric peaks but also extend to radii be- 
yond 2". Although NGC 4889 could host a stellar torus 
with a double-peaked velocity profile, the torus would 
only contribute ~^ 1% of the total signal in the corre- 
sponding LOSVDs, too little to be fully responsible for 
the LOSVDs' non-Gaussianity. 
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Fig. 5. — Stellar kinematics in NGC 4889. Panels (a), (c), (e), and (g) are 2-dimensional maps from GMOS IFS data. The horizontal 
dashed line represents the photometric major axis. Panels (b), (d), (f), and (h) show kinematic moments as a function of radius, after 
averaging kinematic mo ments from d ifferent polar angles in each quadrant of NGC 4889. Panels (b) and (d) also include major-axis 
kinematic moments from lLoubser et~al.. (,2008 '). For the 1-dimensional plots, the values of D^ad and h^ have been inverted on the west side 
of the galaxy, (a) and (b): radial velocity, (c) and (d): line-of-sight velocity dispersion, (e) and (f): Gauss-Hermite moment /13. (g) and 
(h): Gauss-Hermite moment h^. 



The double-peaked feature in our residual image of 
NGC 4889 appears to be qualitatively similar to the cen- 
tral structure seen in a number of high-luminosity galax- 
ies, such as NGC 4073, NG C 6876, and the BCGs N GC 
910, IC 1733, and IC 4329 (jLauer et alJ 120021 12005[) . In 
these galaxies the surface brigh tness profile actual ly has 
a local minimum at the center. iLauer et al.l (|20Q2f) have 
argued that such an appearance could be due to the ad- 



dition of a diffuse torus to the cores of the galaxies, or 
could alternatively be evidence for the evacuation of stars 
from the center by a merging pair of supermassive black 
holes bound in a tight binary. In the case of NCG 4889, 
the torus is not strong enough to create an apparent sur- 
face brightness minimum at the center of the galaxy, and 
would thus be a weaker example of the phenomenon seen 
in the ILauer et al.l (|2002l I2005D galaxies. The radius of 



the torus in NGC 4889 is ~ rinf- In simulations by 
iZier fc Biermaniil (|2001f) . a black hole binary with ini- 
tial separation r-mf scours a torus of radius ~ 3 rinf on a 
timescale comparable to the hardening and coalescence 
of the binary. 

Alternative scenarios for asymmetric kinematics in- 
clude an ongoing minor merger, a surviving nucleus from 
a cannibalized satellite galaxy, or a chance alignment 
with a foreground satellite. However, none of these sce- 
narios are consistent with the highly regular photomet- 
ric contours in NGC 4889: even the 1% deviations from 
perfect ellipses are distributed symmetrically about the 
galaxy center. Another post-merger scenario is a black 
hole displaced from the center by gravitational wave re- 
coil. Simulations of merger remnants with recoiling black 
holes suggest that even with a modest recoil velocity, a 
black hole can maintai n an observable displacement in 
a gas - poor galaxy fe.g..plecha et aLll2011t IGuedes et al.l 



1201 It iSiiacki et al.ll201lD . and damping could be partic- 
ularly weak with a shallow stellar mass profile like the 
core of NGC 4889. However, if a recoiling black hole af- 
fected a sufficient number of stars to displace the galaxy's 
velocity dispersion peak, we would expect to observe a 
similarly prominent photometric disturbance, or even a 
cusp of stars surrounding the black hole. 

Although we do not find an obvious explanation for our 
asymmetric velocity dispersion measurements in NGC 
4889, we suggest that they arise from an asymmetric 
orbital structure within a largely symmetric spatial dis- 
tribution of stars. In Section 15.31 ^^ illustrate how tan- 
gential orbital bias in the vicinity of the black hole can 
produce a central minimum in the line-of-sight velocity 
dispersion. A stellar torus is consistent with excess tan- 
gential orbits as well as our photometric data. However, 
a circular or slightly eccentric torus of stars cannot fully 
explain the velocity dispersion features in the southeast 
and northwest quadrants of NGC 4889 (see Fig. [St). 
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Fig. 6. — LOSVD near the velocity dispersion peak in NGC 4889, 
corresponding to a radius of 1.4" (700 pc) along the east side of 
the major axis. The solid black line is the non-parametric LOSVD, 
and the dashed blue line is the best-fit Gaussian profile. The thick 
red line represents the residual velocity profile after the Gaussian 
profile is subtracted. The residual profile includes tangential or- 
bits centered near ±450kms~^. Negative features in the residual 
profile represent an under-abundance of stars at the corresponding 
line-of-sight velocities. 



4.2. NGC 3842 

Figure [S] shows two-dimensional maps of the kinematic 
moments in NGC 3842, as measured by CMOS, as well 
as radial profiles of each moment as measured by CMOS 
and the Mitcheh Spectrograph. NGC 3842 exhibits the 
simplest kinematics of the four BCGs. There is negligible 
rotation, and a increases virtually monotonically toward 
the center. The maps of /13 and /14 are largely feature- 
less, and /i4 is slightly negative on average. The four 
quadrants of NGC 3842 are similar enough to average 
their kinematics and model a single set of LOSVDs. The 
kinematics from the Mitchell Spectrograph are consistent 
with a continuation of the trends observed at small radii 
with CMOS. 

4.3. NGC 7768 

Figure [9] shows two-dimensional maps of the kinematic 
moments from OSIRIS observations of NGC 7768, and 
radial profiles of each moment from OSIRIS and the 
Mitchell Spectrograph. Unlike most BCGs, NGC 7768 
exhibits a strong rotation curve along the major axis, 
with |wrad| ^ 100 kms^^ at radii beyond 0.2". We mea- 
sure a = 322 ±32 km s~^ in our central OSIRIS bin, out- 
side of which a drops suddenly and remains near 230- 
260kms~^ over the full extent of our radial coverage. 
We measure h^ to be anti-correlated with Wradj particu- 
larly for OSIRIS data. This anti-co rrelation is common 
in rotating early-type galaxies (e.g.. iBender et ai]|1994l : 
IKrainovicet"ari[201ll) . OSIRIS and the Mitcheh Spec- 
tr ograph both measure boxy LOSVDs (/14 < 0). 

ILoubser et all ^(M\ observed NGC 7768 with ISIS 
on WHT, measuring major-axis kinematics out to 11.4". 
At radii of 0.1-1", we measure a steeper radial veloc- 
ity gradient and lower velocity dispersion values with 
OSIRIS. The differences between our measurements and 
those of Loiibscrct aL (2008) are reconciled by consid- 
ering the respective spatial resolutions of the OSIRIS 
and WHT obser vations: whereas O SIRIS is nearly 
diffraction-hmited. ILoubser et al.l (|2008[ ) reported 1" see- 
ing at WHT. Although the poorly resolved velocity gra- 
dient artificially increases the WHT velocity dispersion 
measurements, the quadratic sum v'^^^ + a^ is consistent 
for OSIRIS and WHT data. Only OSIRIS detects the 
central rise to cr > 300kms~^. 

Similarly, we attribute the low value of v^ad in our 
innermost Mitchell Spectrograph bin (Figure IHb) to an 
overlap between the corresponding fiber and the central 
velocity dispersion gradient. We have included this data 
point in stellar orbit mo dels of NGC 7768, a, long with 
five measurements from ILoubser et al.l (|2008f ) at radii 
between 1" and 3". This combination of data is illus- 
trated in the radial plots in Figure [HI In a separate trial, 
we replaced the innermost LOSVD from the Mitchell 
Spectrograph with t hree additional measurements from 
ILoubser et all (|2008n . covering radn of 3.0-5.1". All trials 
yielded consistent values of M, and Af^/L^;. 

4.4. NGC 2832 

Figure [10] shows two-dimensional maps of the kine- 
matic moments from CMOS observations of NGC 2832, 
as well as radial profiles from CMOS and the Mitchell 
Spectrograph. Unlike NGC 4889 and NGC 3842, which 
exhibit either a central or off-center local maximum in 
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Fig. 7.— (a) Residual image of NGC 4889, after fitting elliptical isophotes to a de-convolved image from HST/WFPC2 (F606W). (b) 
Same image, with contrast adjusted to highlight the ~ 1% excess light from a possible stellar torus. The dashed white line indicates the 
major axis of NGC 4889. The white circles are each centered on the major axis 1.4" from the galaxy center. Each circle has a radius of 
1.2". The compact bright spots in both panels are globular clusters. 



a, our velocity dispersion map of NGC 2832 is nearly 
featureless, with values near 360 km s"^ out to r sa 3". 
CMOS IFS data for NGC 2832 show some evidence of 
rotation from north to south, although this signal is 
strangely absent on the east side of the galaxy. There 
is a similar discrepancy in our measurements of h^: we 
measure /13 < on the east side and /13 > on the west 
side. 

ILoubser et aP E(M\ observed NGC 2832 with CMOS 
in long-slit mode, with the slit oriented 67° from the 
galaxy's major axis. They detected a kinematically de- 
coupled core, characterized by rotation out to approx- 
imately 4" along the slit axis. Our measurements of 
^rad agree with the long-slit data on the west side of 
the galaxy, but do not follow the same trend on the east 
side. IFS and long-slit measurements both detect a flat 
velocity dispersion profile for r < 3", with similar aver- 
age values of a. We measure significantly lower velocity 
dispersions with the Mitchell Spectrograph, in the range 
-260-320kms~^ 

The majority of our Mitchell Spectrograph measure- 
ments are at r > 7" and do not overlap with the 
CMOS IFS fie l d of v iew. We use the long-slit data from 
ILoubser et alj ()2008D to help bridge the apparent gap be- 
tween the CMOS and Mitcheh Spectrograph IFS data. 
For radii between 4" and 7" we have averaged multiple 
long-slit data points from ILoubser et aTl (|2008D and pro- 
duced three LOSVDs, which we include in our stellar 
orbit models of NGC 2832. The corresponding values of 
Urad and (7 are included in Figures [TOb and llOH . 

Although measurements by Loubser et al. (2008) indi- 
cate that a may indeed fall from ^ 360knis~^ near the 
center of NGC 2832 to ~ 280 km s^^ at r ~ 10", addi- 
tional discrepancies point to possible systematic errors 
in our kinematics from CMOS and/or the Mitchell Spec- 
trograph. Our LOSVDs from CMOS are peaky (h^ > 0), 
whereas our LOSVDs from the Mitchell Spectrograph are 
boxy (/14 < 0). Furthermore, the two Mitchell Spectro- 
graph fibers nearest the galaxy center yield significantly 



lower velocity dispersions than CMOS IFS or long-slit 
data at overlapping radii. For all galaxies, we exclude the 
innermost Mitchell Spectrograph data from our models, 
so as not to dilute the more highly resolved CMOS data. 
Stellar template mismatch or errant treatment of the 
galaxy and stellar template spectra can bias b oth a 
and h d, (e.g.,[Rix fc White 1992; van der Marcl & Franxl 
[T993I: [Carter elT all ITOM lEmsellem et all 12001 . Al- 
though we have attempted to treat CMOS IFS and 
Mitchell Spectrograph data consistently for each of the 
three BCGs, NGC 2832 exhibits a particularly large 
number of internal inconsistencies. This is considered 
in our final assessment of stellar orbit models for NGC 
2832 (Section EXl]). 

5. STELLAR ORBIT MODELS AND BLACK HOLE MASSES 
5.1. Axisymmetric Orbit Models 

We generate stellar orbit models of eac h galaxy using 
the st atic potential method introduced by iSchwarzschildl 
(|1979D . We use the axisymmetric modeling algorithm 
described in detail in Gebhardt et al. (2000b ; 2003), 
Thomas et al. (2004; 2005), ai id iSiopis et all (12009.) 
Similar models ar e prese nt ed in iRichstone fc Tremaing 
(I1984D. iRix et al.l ([l99 7l) Ivan der Ma rel et al.l (I1998D. 
iCretton et al.l (|1999[ ). and lValluri et a l. (2004). 

We assume that each galaxy includes three mass com- 
ponents - stars, a central black hole, and an extended 
dark matter halo - described by the radial density pro- 
file 

P{r) = ^Hr) + M.5{r) + Phaio(r) , (3) 

where i>{r) is the observed luminosity density (see Fig- 
ure [1]) . The resulting gravitational potential depends on 
M,, the stellar mass-to-light ratio M^/L {AI^^/Ln for R- 
band photometry), and two parameters describing the 
dark matter density profile /C»haio(7")- Our models are sym- 
metric about the z-axis (corresponding to the projected 
minor axis) and the equatorial plane (z = 0). Stellar or- 
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Fig. 8. — Stellar kinematics in NGC 3842. Panels (a), (c), (c), and (g) are 2-dimensional maps from GMOS IFS data. The horizontal 
dashed line represents the photometric major axis. Panels (b), (d), (f), and (h) are radial plots, including GMOS and Mitchell Spectrograph 
data. In order to construct the radial plots, the kinematic moments from different quadrants and polar angles have been averaged, after 
inverting the values of tij-ad and /13 in the southeast and southwrest quadrants, (a) and (b); radial velocity, (c) and (d): line-of-sight velocity 
dispersion, (e) and (f): Gauss- Hermite moment ^3. (g) and (h): Gauss- Hcrmite moment /14. 



bits are generated by propagating test particles through 
the potential and computing their time-averaged veloci- 
ties in each bin in a polar grid. Typical models produce 
~ 30, 000 bound orbits, including both signs of the an- 
gular momentum component Lz- 

To match each observed LOSVD, we compute a model 
LOSVD from the projected velocities of individual orbits 
in the spatial region corresponding to the extracted spec- 
trum and assumed PSF. Each orbit is assigned a scalar 



weight, and the orbital weights are varied to optimize 
the fit between the observed and model LOSVDs. The 
goodness of fit for each model is characterized by the x^ 
parameter: 



N, 
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where Ci data and Ci model are LOSVDs in each of the 
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Fig. 9. — Stellar kinematics in NGC 7768. Panels (a), (c), (e), and (g) are 2-diniensional maps from OSIRIS IFS data. The horizontal 
dashed line represents the pho tom etric major axis. Panels (b), (d), (f), and (h) are radial plots, including OSIRIS and Mitchell Spectrograph 
data, plus long-slit data from [Loubscr ct al. (2008) at radii from 1" to 3". In order to construct the radial plots, the kinematic moments 
from different hemispheres and polar angles have been averaged, after inverting the values of ii^ad a^nd h-3 on the west side of the galaxy, (a) 
and (b): radial velocity, (c) and (d): line-of-sight velocity dispersion, (e) and (f): Gauss-Hermite moment hg. (g) and (h): Gauss-Hermite 
moment h4. 



i = 1, ...Nb spatial bins, and ^^{vj) is the squared un- 
certainty in £i_data at velocity bin Vj. The weights are 
constrained such that the summed spatial distribution of 
all weighted orbits must match the observed luminosity 
density profile. 

For each galaxy, we have run multiple trials of orbit 
models. Each trial assumes a particular dark matter halo 
profile and fits a specific combination of kinematic data 
to several hundred models covering a finely sampled grid 



in M, and M^^/L. The resulting confidence limits in M, 
and Mi^/ L are computed by analyzing the distribution of 
X^ val ues, according to the cumulative likelihood method 
of McConnell et al.l (|2011aD . Our modeling trials for the 
four galaxies are summarized in Tables |3]|6l 

High-resolution and wide-field data play complemen- 
tary roles in constraining the gravitational potential of 
each BCG. Kinematic measurements from CMOS are 
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Fig. 10. — Stellar kinematics in NGC 2832. Panels (a), (c), (e), and (g) are 2-dimensional maps from GMOS IFS data. The horizontal 
dashed line represents the pho tometric major axis. Panels (b), (d), (f), and (h) are radial plots, including GMOS and Mitchell Spectrograph 
data, plus long-slit data fromXpubscr ct al. (200^) at radii from 4" to 7". In order to construct the radial plots, the kinematic moments 
from different quadrants and polar angles have been averaged, after inverting the values of fj-ad and /13 in the southeast and southwest 
quadrants, (a) and (b): radial velocity, (c) and (d): line-of-sight velocity dispersion, (e) and (f): Gauss-Hermite moment h-j. (g) and (h): 
Gauss-Hermite moment ^4. 



sufficient to detect a black hole in NGC 4889 and NGC 
3842, in part because of excellent seeing. Still, the 
black hole's gravitational signature is not fully decou- 
pled from the enclosed stellar mass, and improved mea- 
surements of Mi,/L yield tighter confidence intervals in 
M,. Wide-field kinematics from the Mitchell Spectro- 
graph and long-slit instruments constrain the relative 
contributions of stars and dark matter and increase the 



accuracy and precision of Mi,/ L and M,. 

In principle, the orbit modeling code can accommo- 
date any dark matter density profile. We have restricted 
our experiments to two f unctional forms: the core d LOG 
profile (descri bed in, e.g..lMcConnell et al.ll2011aj) . or the 
NFW profile (jNavarro et al.lll99(][) . The free parameters 
in the LOG profile are the asymptotic circular speed Vc 
and the core radius Vc- The free parameters in the NFW 



13 



profile ar e the concentration p arameter c and the scale 
radius r^ . iThomas et al.l (|2007[ ) used long-slit kinematics 
to determine the best-fitting LOG and NFW profiles for 
NGC 4889. Our models of NGC 4889 approximate their 
best-fitting LOG profile. 

The literature does not provide a robust estimate of the 
dark matter halo profile for NGC 3842, and our kine- 
matics from the Mitchell Spectrograph extend to only 
O.Qreff. We have tested LOG profiles with a few difi^er- 
ent values of Vc, as the total halo mass is proportional to 
w^. After determining the best value of Vc (350 km s^^, 
with Tc = 8.0 kpc), we tested a comparable NFW halo, 
scaled to enclose the same mass within the outer radius 
of our kinematic data. We required c and r^ to follow 
the relation 



3 X 1013 M, 







200f pcrit C3 



10oT5(l-''^~'°Sl"'^) 



(5) 



(|Rix et al.lll997t iMcConneh et al.ll2011al) . where pcrit is 
the present-day critical density. The resulting NFW halo 
has c = 13.5 and rg = 31.2 kpc, and fits our data to the 
same confidence level as the LOG profile (Axmin = 0-4) ■ 
Even though the density of the NFW profile increases 
toward r — 0, the total mass near the center is dominated 
by stars and the supermassive black hole. We further 
describe our results modeling different dark matter halos 
in Section [5X2] 

Our Mitchell Spectrograph kinematics extend to 
1.2rcff for NGC 7768 and 1.5rcff for NGC 2832, but in 
both cases fail to tightly constrain the dark matter pro- 
file. In Section 15.2.31 we describe our treatment of dark 
matter in NGC 7768 and note that our measurement 
of M, is insensitive to the dark matter profile, in Sec- 
tion [5231 we describe our dark matter models for NGC 
2832 and how they affect our measurement of M, . 

5.2. Black Hole Masses and Mass-to-Light Ratios 

As in iMcConnell et al.l (|2011af ). we have determined 
confidence intervals in M, and M^,/ Lji by running nu- 
merous models and integrating the relative likelihood 
function, P oc e~2(x -Xmin). Results for individual galax- 
ies are discussed below. 

5.2.L NGC 4889 

In order to address the dramatic asymmetries in the 
stellar kinematics near the center of NGC 4889, we have 
independently run models to fit each observed quadrant 
of the galaxy, assuming the same central LOSVD (r < 
0.25") for all quadrants. For each quadrant, we fit IFS 
kinemati cs from CMOS and long-slit (WHT) kinemat- 
ics from iLoubser et al.l (|2008f ). The resulting contours 
in x^{M,,M^,/Lr) are illustrated in Figure El There 
is a large degree of overlap between the 2-dimensional 
68% confidence regions for the northeast, southeast, and 
northwest quadrants, each of which contain a local max- 
imum in a. After marginalizing over M^,/ Lr, these three 
quadrants yield mutually consistent black hole mass es- 
timates. In the southwest quadrant, a decreases nearly 
monotonically with increasing r, and the marginalized 
68% confidence interval M, = 5.5 - 17.0 x 10^ M© only 
agrees with the northeast quadrant. We wish to repre- 
sent the constraints that all four quadrants impose upon 



the central black hole mass, and so we adopt the most ex- 
treme range of confidence limits, M, = 5.5 — 37x 10^ M©. 
In other words, we exclude only those models whose or- 
bit solution is an outlier in all four quadrants. We define 
the best-fit black hole mass as the midpoint of the confi- 
dence interval above, such that A/, « 21 ± 15 x 10^ M©. 
The stellar mass-to-light ratios from all four quadrants 
agree within 68% confidence limits. The extreme 68% 
confidence interval is Mi,/Lf( = 4.2 — T.QMqLZ^^. 

For a torus of stars in NGC 4889, we can estimate 
the total enclosed mass as Mcnc = rv^ jG. We consider 
a characteristic velocity, f, of 450 km s^^, based on the 
double-peaked residual features from LOSVDs near the 
major axis, and an average radius, r, of 1.4" (700 pc). 
These estimates yield Mc„c = 3.3 x lO^^Mg. Our stel- 
lar orbit models of NGC 4889 find M^jL^. in the range 
of 4.2 — 7.6 MqLq^^, corresponding to stellar masses of 

0.7- 1.2 X IO^^Mq within 1.4". The remaining mass of 
2.1 — 2.6 X 10^° Mq lies near the middle of our adopted 
confidence interval for M,. 

Although NGC 4889 does not show strong evidence for 
an off-center black hole, we wish to estimate the range 
of black hole masses that could reside at the location 
of the global maximum in u. We have fit models to 
the LOSVDs from the east side of NGC 4889, spatially 
offset such that the center of the axisymmetric gravita- 
tional potential is aligned with the LOSVD at r = 1.4" 
along the major axis. These models produce a 68% con- 
fidence interval M, = 8.7 - 24 x lO'' Af©, which falls 
entirely within the range spanned by the models of indi- 
vidual quadrants. Although these "recentered" models 
do not accurately reflect the spatial relationship between 
the stellar mass distribution and the observed kinemat- 
ics, they support our overall conclusion that a black hole 
of at least 5.5 x 10^ Mq is required to produce the ob- 
served line-of-sight velocity dispersions near 400 km s"^. 
Although the observed line-of-sight velocity dispersion 
decreases toward the center of NGC 4889, the best-fitting 
stellar orbit model for each quadrant predicts that the 
3-dimensional velocity dispersion ctsd rises toward the 
black hole. A non-isotropic orbital distribution, biased 
toward tangential orbits, causes the projected velocity 
dispersion to decrease even as ctsd increases. This is fur- 
ther discussed in Section [Ol and illustrated in Figure [T71 

5.2.2. NGC 3842 

Models using our best estimate of the dark matter halo 
in NGC 3842 yield 68% confidence intervals of M, = 
7.2-12.7x10^ Mq and M^/ Lr = AA-5.8 MqL-^j^. The 
confidence interval for M, overlaps with the confidence 
intervals from all other trials, even models with no dark 
matter. The best-fit value of AT, decreases by 26% when 
dark matter is omitted and increases by 46% for our most 
massive halo. We observe a similar level of sensitivity 
in Mi,/Lji, where the most extreme trials yield best-fit 
values of 3.7 and 7.2AfQLQ"'^. Figure [T2| illustrates the 

X^ contours for models with our best-fitting NFW profile 
and with no dark matter. 

Before acquiring Mitchell Spectrograph data, we mod- 
eled NGC 3842 with CMOS IFS dat a plus major-axis 
kinematics from 'Loubs er et al.l (|2008D , measured using 
ISIS at WHT. The more recent Mitchell Spectrograph 
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TABLE 3 

AXISYMMETRIC MODELS OF NGC 4889 



Data 



(km s -"-) (kpc) 



Afh{20 kpc) 
(10" Mq) 



fDMiXeff) 



M. 

(109 






A'dof 



(1) 



(2) 



(3) 



(4) 



(5) 



(6) 



(7) 



M W 



CMOS southwest + WHT 
CMOS southeast + WHT 
CMOS northwest + WHT 
CMOS northeast + WHT 

CMOS west + WHT 
CMOS east + WHT 



425 


8.0 


7.21 


0.54 


425 


8.0 


7.21 


0.57 


425 


8.0 


7.21 


0.56 


425 


8.0 


7.21 


0.55 


425 


8.0 


7.21 


0.54 


425 


8.0 


7.21 


0.58 



9.8 



+7.2 



2&+_l 

27tf 

7+8 



n-t 



12 
29 



+8 

5.5 
+5 



6.6 
5.6 
5.8 
6.1 

6.4 
5.4 



+ 1.0 
1.3 
1.1 
1.4 
1.2 
1.4 
1.2 
1.5 

+ 1.0 
1.2 

+ 1.0 
0.9 



176.6 345 
137.0 420 

170.7 345 
119.4 420 

180.4 330 

217.4 442 



CMOS recentered + WHT 



425 



8.0 



7.21 



0.54 



15-* 



6.5Z 



239.5 442 



Notes: Column 1: Data sets included in trial. Column 2: circular velocity of LOG dark m atter halo. Column 3 : core radius 
of LOG dark matter halo. For NGC 4889, we have approximated the best-fitting halo from I Thomas et al.l (120071 ). Column 4: 
enclosed halo mass within 20 kpc. The "no dark matter" case has v^ = 0.01 km s~^, r^ = 20.0 kpc, and Mh ~ 200 M©. Column 
5: dark matter fraction within one effective radius, assuming the best-fit values of M, and Mi^ILn- Column 6: best-fit black 
hole mass. Errors correspond to 68% confidence intervals. Column 7: best-fit i?-band stellar mass-to-light ratio. Quoted errors 
correspond to 68% confidence intervals. Column 8: minimum ^ value for all models. Column 9; degrees of freedom in model 
fits to LOSVDs. Computed values include a smoothing factor of 1 degree of freedom per 2 velocity bins for non-parametric 
LOSVDs from CMOS. 
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Fig. 11.— Contours of x^ versus kh/LR and M., for models of NGC 4889, fitting LOSVDs from the CMOS IFS and fromToub ser et al.l 
II2008I ). Data from each quadrant of the galaxy were fit independently, and the respective results are shown in panels (a) through (d). 
Contours of Ax^ = 1, 2.71, and 6.63 represent confidence levels of 68%, 90%, and 99% for one free parameter. Small crosses denote 
individual models. 



data yield slightly larger values of M,, but are con- 
sistent within 68% confidence limits. We prefer using 
the Mitchell Spectrograph data because they extend to 



larger radii (35.3" versus 20.8") and provide full two- 
dimensional spatial sampling. 
OSIRIS and CMOS provide independent measure- 
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TABLE 4 

AXISYMMETRIC MODELS OF NGC 3842 



Data 



(1) 



(km s -*-) (kpc) 



(2) 



(kpc) 



(3) (4) (5) 



Mh(20 kpc) 
(10" Mo) 

(6) 



/DM(»'cff) 



(7) 



M. 

(109 

(8) 



M./Lr xl 



(9) 



A^do 



(10) (11) 



CMOS + Mitchell 
CMOS + Mitchell 
CMOS + Mitchell 
CMOS + Mitchell 

CMOS + OSIRIS + Mitchell 
CMOS + OSIRIS + Mitchell 



CMOS 
CMOS 



WHT 
WHT 



0.01 


20.0 


2 X 10"^ 





7 9+2.1 
'•^-2.9 


350 


8.0 


4.89 


0.51 


10.2+1? 




13.5 31.2 


5.09 


0.53 


9.7tli 


500 


8.0 


9.98 


0.75 


13.2t|J 


0.01 


20.0 


2 X 10"^ 





6^3l}1 


500 


8.0 


9.98 


0.74 


10.2tlJ 


0.01 


20.0 


2 X 10^'J 





6.8+- 


500 


8.0 


9.98 


0.71 


9.9l^- 



7.2 
5.4 
5.1 
3.7: 

7.1 
4.0 

7.3 

4.8 



+0.6 
0.5 


92.6 


345 


+0.6 
-0.7 
+0.7 

+8:1 

-0.6 


89.3 
88.9 
96.2 


345 
345 

345 


+0.6 

+0:4 
-0.5 


129.6 
138.4 


368 
368 


+0.5 
-0.6 
+0.5 
-0.8 


355.6 
349.3 


525 
525 



Notes: Column 1: Data sets included in trial. Column 2: circular velocity of LOG dark matter halo. Column 3: core radius of 
LOG dark matter halo. Column 4; concentration parameter for NFW dark matter halo. Column 5: scale radius for NFW dark 
matter halo. Column 6: enclosed halo mass within 20 kpc. Our best-fitting halo is an NFW profile with c = 13.5 and Va = 3L2 



kpc. The equivalent LOG halo {vc = 350 km s 

-1 



— 8.0 kpc) fits our data to the same level of confidence (Axmin = 0.4). The 
"no dark matter" case has Vc — 0.01 km s~'^, r^ = 20.0 kpc, and Mh ~ 200 Mq. Omitting dark matter from the models causes 
M, to decrease by 26%, consistent within errors. Column 7: dark matter fraction within one effective radius, assuming the 
best-fit values of M, and Mi,/LR. Column 8: best-fit black hole mass. Errors correspond to 68% confidence intervals. Column 
9: best-fit 7?-band stellar mass-to-light ratio. Quoted errors correspond to 68% confidence intervals. Column 10: minimum x^ 
value for all models. Column 11; degrees of freedom in model fits to LOSVDs. Computed values include a smoothing factor of 
1 degree of freedom per 2 velocity bins for non-parametric LOSVDs from GMOS and OSIRIS. 
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Fig. 12.— Contours of x^ versus M*/Lfl and M., for models of NGC 3842, fitting LOSVDs from GMOS and the Mitchell Spectrograph. 
(a) Results for our best-fitting NFW dark matter profile (c = 13.5, rs = 31.2 kpc). (b) Results for models without dark matter. Contours of 
Ax^ = 1, 2.71, and 6.63 represent confidence levels of 68%, 90%, and 99% for one free parameter. Small crosses denote individual models. 



ments of stellar kinematics in NGC 3842, for r < 0.7". 
In spite of the potential for higher spatial resolution, sys- 
tematic contaminants in our OSIRIS spectra force us to 
increase S/N by binning the data to identical spatial 
scales as for GMOS. We have run orbit models fitting 
LOSVDs from OSIRIS and GMOS simuhaneously (as 
well as Mitchell Spectrograph data at large radii). In- 
cluding the OSIRIS data causes the best-fit value of M, 
to decrease by up to 23%, likely because OSIRIS data 
show a less drastic increase in a than data from GMOS. 
For the best-fit models, including OSIRIS data produces 
higher average x^ values per LOSVD, even after ignor- 
ing the central regions where LOSVDs from OSIRIS and 
GMOS arc not fully consistent. The kinematic fits to 



OSIRIS data may have systematic errors, resulting from 
imperfect sky subtraction and difficulty defining the H- 
band continuum. Consequently, we judge the models of 
NGC 3842 with only GMOS and MitcheU Spectrograph 
data to be the most reliable. 

5.2.3. NGC 7768 

For NGC 7768, we have combined IFU data from 
OSIRIS and the Mitchell Spectrograph, plus five long-slit 
measurements at radii of 1-3". We have run two trials 
with full sampling of M, and Mi,/L. One trial includes a 
LOG halo profile, matching our best-fitting approxima- 

= 8.0 kpc). The 
The black hole 



-1-. 



tion for NGC 3842 [v^ = 350 km s' 
other trial docs not include dark matter. 
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TABLE 5 

AXISYMMETRIC MODELS OF NGC 7768 



Data 



(1) 



(km s -'-) (kpc) 



(2) 



(3) 



Mh(20 kpc) 
(10" Mq) 

(4) 



/DM(»'cff) 



(5) 



M. 
M 



M^/Lr 
(7) 



A^do 



M W 



OSIRIS + Mitchell + WHT 0.01 

OSIRIS + Mitchell + WHT 350 



20.0 2 X 10"^ 
8.0 4.89 




0.32 



1.3 
1.2 



+0.5 
0.4 

+0.5 
0.4 



6.1 
5.2 



+0.7 
0.8 

+0.5 
0.6 



77.4 
78.8 



285 
285 



Notes: See Table [3] for notes and definitions. For WHT measurements, we only include five LOSVDs, covering radii between 
1" and 3". The LOSVDs extracted from the Mitchell Spectrograph cover five radial bins between 3.0" and 28.5" (0.1 - 1.2 
X reft). Several models with different dark matter profiles, including no dark matter and the dark matter model listed here, all 
have Ax^ < 1- Our best-fit black hole mass for NGC 7768 is insensitive to the assumed dark matter profile, a testament to the 



strength of our OSIRIS data. 
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Fig. 13. — Contours of x^ versus Ivh/La and M., for models of NGC 7768, fitting LOSVDs from OSIRIS and the Mitchell Spectrograph, 
plus five long-slit data points from lLoubser et al.l II2008I) . (a) Results for a LOG dark matter profile {vc = 350 kms~^, re = 8.0 kpc). (b) 
Results for models without dark matter. Contours of Ax^ = 1, 2.71, and 6.63 represent confidence levels of 68%, 90%, and 99% for one 
free parameter. Small crosses denote individual models. 



masses from the two trials are in excellent agreement, 
and the best-fit values of Mi^/Ln agree within 68% con- 
fidence limits. The trial with no dark matter actually 
yields a lower value of Xmin' albeit with marginal sig- 
nificance (Ax^jjj = 1.4). Considering the lower value 
of Xmin from the trial without dark matter, we formally 
adopt a black hole mass of 1.3lo!4 x 10^ M©. The range 
of M^/Lr from both trials is 4.6-6.8 MqL~\ (68% con- 
fidence). The respective x^ contours are shown in Fig- 
ure [131 

Given our initial trials' consistent measurements of A/, 
and Mi,/ L^ we have not run finely sampled trials for ad- 
ditional dark matter profiles. In order to check whether 
any dark matter profile for NGC 7768 can fit our kine- 
matics as well as the trial without dark matter, we have 
run a grid of models coarsely sampling M^/L, Vc, and 
Tc- This trial yields several models whose x^ values are 
consistent with Xmin for the trial without dark matter. 



5.2.4. NGC 2832 

For NGC 2832, we have combined IFU data from 
CMOS and the Mitchell Spectrograph, along with three 
long-slit measurements at radii of 4-7", and have run 
full trials varying M, and M*/L for two dark matter 



halos. Figure [U illustrates the x^ contours, and Fig- 
ure [T5] shows the likelihood distribution of M, from each 
trial, after marginalizing with respect to Mi,/L. Models 
for the less massive halo {vc = 350 km s^^; re = 32.0 
kpc) yield a best-fit black hole mass of 2.9 x 10^ M© and 
a 68% upper confidence limit of 5.2 x 10^ M©. How- 
ever, the likelihood distribution falls slowly as Af, ap- 
proaches zero, and has sufficient noise between M, = 
and A'f, = 6.0 x 10^ Mq to cast doubt upon any lower 
limit in M,. In contrast, the more massive halo {vc = 560 
km s~^; Tc = 40.0 kpc) yields a best-fit black hole mass 
of 5.9 x 10^ Mq, and the corresponding likelihood distri- 
bution declines cleanly to both sides, producing a 68% 
confidence interval of M, = 2.7-9.0x10^ Af©. We obtain 
consistent results for this halo when the tliree long-slit 
data points at intermediate radii are excluded. 

When data from CMOS and the Mitchell Spectrograph 
are modeled simultaneously, the minimum x^ values for 
all input Af, and M^^/L significantly favor the more mas- 
sive halo (Ax^jjj = 3.9). However, there are multiple 
reasons to be skeptical of the corresponding confidence 
limits for A/,. First, the sharp jump between the velocity 
dispersions measured with CMOS and those measured 
with the Mitchell Spectrograph suggests that systematic 
errors in the kinematic fitting could be responsible for 
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TABLE 6 

AXISYMMETRIC MODELS OF NGC 2832 



Data 



(1) 



vc re A/h(20 kpc) /DM(reff) M. M^jLn X^i^ 

(kms-i) (kpc) (10" Afe) (10^ M©) (MqL-i^) 



(2) 



(3) 



(4) 



(5) 



(6) 



(7) 



A^H, 



m w 



CMOS IFU + Mitchell + CMOS slit 350 32.0 1.59 0.28 < 5.2 

CMOS IFU + Mitchell + CMOS slit 560 40.0 2.90 0.47 5.9 

CMOS IFU + Mitchell 560 40.0 2.90 0.47 7.2J 



+3.1 
3.2 



8.1 
7.3 



+0.3 
0.5 

+0.4 
0.4 



7.2 



+0.3 



129.0 495 

125.1 495 

115.2 450 



Notes: See Table [3] for notes and definitions. For GMOS slit measurements, we only include three LOSVDs, covering radii 
between 4" and 7". 
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Fig. 14.— Contours of x^ versus Mi^jL n and M., for models of NGC 2832, fitting LOSVDs from GMOS and Mitchell Spectrograph IFS 
data, plus three long-slit data points from lLoubser et al.l {2003). (a-) Results for a high-mass LOG dark matter profile («c = 560 km s~^ , re = 



40.0 kpc). (b) Results for a lower-mass LOG profile (v^ = 350 kms ^^Tc = 32.0 kpc). Contours of Ax^ 
confidence levels of 68%, 90%, and 99% for one free parameter. Small crosses denote individual models. 
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Fig. 15.— Likelihood, P, versus M. for models of NGC 2832, 
after marginalizing with respect to M^/Ljj. The black line rep- 
resents the trial with a high-mass LOG dark matter profile [vc = 
560kms~^,rc = 40.0 kpc), and the grey shaded region represents 
the corresponding 68% confidence interval. The blue line and cyan 
shaded region represent the trial with a lower-mass LOG profile 
(vc = 350kms~^,rc = 32.0 kpc). The diagonally hashed region is 
where the two confidence intervals overlap. Both likelihood distri- 
butions have been normalized to a maximum value of 1. Consid- 
ering both trials, and discrepancies between different instruments' 
kinematic measurements for NGC 2832, we adopt an upper limit 
of M. =9 X 10^ Mq. 



the apparent signature of a massive black hole. Second, 
the best-fit stellar mass-to-light ratio of 7.3 MqL~^^ ex- 
ceeds the upper limits for M^^/Lr in NGC 4889 andNGC 
3842. Third, the best-fit value of M, in NGC 2832 is sen- 
sitive to the dark matter halo profile, and one trial with 
a plausible halo fails to produce a significant black hole 
detection. Given the possibility of systematic errors in 
the kinematics, we do not find convincing evidence for a 
specific dark matter profile in NGC 2832. We conserva- 
tively interpret our models of NGC 2832 to provide an 
upper limit of M, = 9.0 x 10^ Mq, to 68% confidence. 

5.3. Tangentially Biased Orbits 

In NGC 4889, NGC 3842, and NGC 7768, the super- 
massive black hole coincides with a deficiency of radial 
orbits near the galaxy center. In Figure I16[ we plot or- 
bital anisotropy versus radius for each BCG and com- 
pare predictions from our orbital models with the best-fit 
black hole mass and without a black hole. Specifically, 
we examine the ratio of the radial and tangential velocity 
dispersions in the 3-dimensional models. The tangential 
dispersion term includes azimuthal rotation and is de- 
fined as al^^ = i(cr^ + cr| -h W^), such that CTrad/c^tan = 1 

for an isotropic orbit distribution. 

Figure [16] illustrates that for NGC 4889, NGC 3842, 
and NGC 7768, the models with black holes predict 
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strong tangential velocity biases toward the galactic cen- 
ters (solid curves), while the models without black holes 
predict fairly isotropic velocity distributions with little 
radial dependence (dashed curves). After averaging over 
the four quadrants using the best-fit black hole masses, 
the minimum value of cTiad/ctan in NGC 4889 is approx- 
imately 0.3. The best models for NGC 3842 and NGC 
7768 have central values of tTrad/ctan ~ 0.3—0.4. This de- 
gree of tangential bias is comparable to t he central values 
in a handfu l of carly-typc galaxies fc.g.. 'Gebhard t et al.l 
[200l [20071: [Shcn & Gebhardt 2010). NGC 4889 and 
NGC 3842, however, are exceptional in terms of the size 
of the central regions dominated by tangential orbits. In 
both galaxies, this region extends to approximately 2", 
or 1 kpc. By comparison, tangential bias extends to only 
50-200 pc throughout a broad sample of early- type galax- 
ies with black hole mass measurements (Gebhardt et al. 
2003; 2007; Nowak et al. 2008; Shen & Gebhardt 2010; 
Gebhardt et al. 2011; c.f., Giihekin et al. 2009b). In 
NGC 7768, we measure a tangential bias for r < 0.5" 
(270 pc) for a black hole mass of 1 x IO^Mq. 

The scarcity of radial orbits at the centers of massive 
galaxies may be a consequence of past core-scouring by 
one or more pairs of in-spiraling black holes. A'^-body 
simulations of merging galaxies and black holes indi- 
cate that stars on radial orbits are more likely to inter- 
act with a binary black hole, resulti ng in the ir ejection 
from the galaxy cen ter (e.g., Quinlan fc Herng uist 1997; 



Milosavlievic & Mcrritt 2001). Supporting this picture, 
Gebhardt et al. (2003 ) ran stellar orbit models of twelve 



galaxies with black holes and noted that the four core- 
profile galaxies had a higher degree of tangential bias 
than the eight power law galaxies. Photometric data of 
NGC 4889 show a stellar core with a radius of 750 pc 
(jLauer et al. 2007a), similar to the the scale on which 
the tangential bias appears. In the sarn ple of 118 core- 
profile galaxies from iLauer et al.l ()2007aD , only eight ex- 
hibit core radii greater of 750 pc or greater; all eight are 
BCGs or cD galaxies of comparable luminosity. In NGC 
3842, the core radius is 310 pc, roughly one third the size 
of the tangentially biased region. The surface brightness 
profile for NGC 7768 turns over at approximately 220 pc, 
similar to the size of the tangentially biased region. 

It is important to note that the decrease in velocity 
dispersions towards the black hole, as illustrated in Fig- 
ure [5] for NGC 4889, occurs for the observed line-of-sight 
velocity dispersions. The 3-dimensional velocity disper- 
sion (T3D is not directly measureable but is predicted by 
our orbital models. Figure [T71 illustrates the radial trends 
in USD for models of each quadrant of NGC 4889. As ex- 
pectd, (T3D rises towards the center of the galaxy when a 
black hole is included in the model. Figures [16] and [17] to- 
gether show that a non-isotropic orbital distribution, bi- 
ased toward tangential orbits, can cause the light-of-sight 
velocity dispersion to decrease even as ctsd increases. 

Our models of NGC 4889 indicate that orbital 
anisotropy is responsible for the local minimum in a at 
the galaxy center. In NGC 3842, the increase in a is less 
pronounced than would be observed for isotropic orbits 
about a lO^^-M© black hole. For both galaxies, models 
with M, = have enough freedom to imitate the ob- 
served kinematics, but require a nearly isotropic orbital 
distribution to do so (Figure [T6ti. [T6bV In NGC 7768, 
a radial bias is required to reproduce the sharp increase 



in a (Figure [TBb ) . Although the LOSVDs generated by 
models with and without black holes are qualitatively 
similar, the models with M, = produce higher values 
of x^ in every trial. 

NGC 2832 has nearly isotropic orbits at all radii, and 
the orbital discrepancy between models with and without 
a black hole is weaker than in NGC 4889 and NGC 3842. 
This is particularly true for models with the less massive 
dark matter halo, where we lack a confident black hole 
detection. The best-fitting model with the more massive 
halo has M, — 7.2 x 10^ Mq, and crrad/ftan declines to 
0.7 over the the inner 0.3" (150 pc). 

5.4. Possible Systematic Errors 

The common practice for reporting errors in stellar dy- 
namical measurements of Af, and Mi,/ L is to report the 
range of input A/, and M-^/ L values that fall within a 
specific confidence interval for a given set of models. Yet 
modeling a galaxy requires explicit assumptions about 
its stellar content and the structure of its gravitational 
potential. Systematic errors from these assumptions will 
likely increase the overall measurement errors for Af, and 
MJL. 

For each BCG, our models assume an edge-on incli- 
nation. This is indirectly supported by the observed 
axis ratios of 0.68 to 0.86, which imply relatively eccen- 
tric intrinsic shapes even for an edge-on system. Models 
with more face-on inclinations could yield higher black 
hole masses. We also model each galaxy as an oblate 
axisymmetric ellipsoid, in contrast to the presence of 
isophotal twists in NGC 3842 and NGC 2832 (Table P 
and circumstanti al evidence that s o me BCGs are pro- 
late or triaxial (iPorter et all Il99lt iRvden et all 119931: 
^Boylan- Kolchin et al.ll200 6n. In an early experiment with 
triaxial models, the best-fit black hole mass for NGC 
3379 increased by a factor of two, whereas the black hole 
mas s for M32 was unchanged from a xisymmetric mod- 
els (jvan den Bosch fc de Zeeuwll2010D . The increase for 
NGC 3379 may have arisen solely from a change in the 
best-fit inclination (R.C.E. van den Bosch, 2011, private 
communication) . 

Massive elliptical galaxies are known to exhibit 
spati a l gradients in metallicity ie.s... iFisher et "all 
19951: iKobavash i fc Ar iimtol [l999l: iMehlert et al.l 120031: 

m 



Brough et al.l 12007: Gr eene et al.l I2012D . and gradients 



in stellar age and a-element enhancement ha ve been re- 
porte d in some cases (Fisher et al. 1995; Co ccato et al.l 
l2OT0h . 

Massive ellipticals could exhibit corresponding gradi- 
ents in Mi,/ L^ yet our orbit models assume that Mi,/L 
is constant. This assumption could bias our estimates 
of the enclosed stellar mass at large or small radii, de- 
pending upon the steepness of the gradient and the radial 
sampling of our kinematic data. Preliminary single-burst 
stellar population modeling of M87 indicates that Mi,/L 
decreases by ~ 50% from r = to r ~ 2roff (Graves & 
Murphy, 2012, in prep). Most of our kinematics for the 
BCGs herein correspond to r < roff, and we estimate 
that a gradient in M^,/ L would bias our measurement of 
the central Af*/L by < 10%. Kinematic and photomet- 
ric data probe the galaxy's total mass-to-light ratio, and 
models with dark matter have more fiexibility to repro- 
duce each galaxy's total mass profile. Trade-offs between 
dark matter and gradients in Mi,/L could mitigate poten- 
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Fig. 16. — Anisotropy of model stellar orbits in each galaxy, as a function of radius. Solid lines represent the best-fitting models for 
each set of LOSVDs, and dashed lines represent the best-fitting models without a black hole. We define (Trad a-^d ctan such that their 
ratio is unity for an isotropic orbit distribution, (a) NGC 4889. Different colors represent the northeast (black; M, = 1.2 X IO'^'^Mq), 



^'^ Mq), and southwest (green; M. = 6.0 X 10^ M©) quadrants. 
X 10^ Mq. (c) NGC 7768. The best-fitting model with dark 



northwest (purple; M, = 2.6 X 10^" Mq), southeast (red; M, = 2.2 X 10 

(b) NGC 3842. The best fitting model has c = 13.5, r^ = 31.2 kpc, M, 

matter (black lines) has Vc = 350kms~^, re = 8.0 kpc, and M. = 1.0 X lO" Mq. The best-fitting model without dark matter (cyan lines) 

has M. = 1.0 X 10^ Mq. (d) NGC 2832. The best-fitting model for the high-mass dark matter halo (black lines) has Vc = 560 kms'^, 

re = 40.0 kpc, and M, = 7.2 X 10^ Mq. The best-fitting model for the low-mass halo (cyan lines) has Vc = 350kms~^, r^ = 32.0 kpc, and 

M. = 8.0 X 10** Mq. 
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rants. 



tial biases in our measurements of A/,. Still, our knowl- 
edge of the individual mass components in BCGs would 
be improved by using stellar line indices and population 



models to measure gradients in Mi,/L independently. 

A fundamental assumption of all orbit superposition 
models is that the stellar motions reflect a steady-state 
gravitational potential. The strongly asymmetric kine- 
matics in NGC 4889 present the possibility of an unre- 
laxed stellar system, although we do not witness corre- 
sponding photometric irregularities. Even with our con- 
servative treatment of the asymmetric data, our models 
would be unable to assess the central black hole mass if 
the kinematics of NGC 4889 reflected a temporary con- 
dition. 

6. CONCLUSION AND DISCUSSION 

We have reported three black hole masses (NGC 4889, 
NGC 3842, NGC 7768) and one upper limit on M. 
(NGC 2832) in four of the nearby Universe's most mas- 
sive galaxies at distances of ^ 100 Mpc. As the ba- 
sis for our analysis, we have presented high-resolution, 
2-dimensional measurements of stellar kinematics using 
CMOS at Gemini Observatory North and OSIRIS at 
Keck Observatory. For three of our galaxies, we have 
used wide-field stellar kinematics from the Mitchell Spec- 
trograph at McDonald Observatory. 

We have derived black hole masses from stellar kine- 
matics using orbit superposition models. We find M, = 
2.1±1.6x IQi" Mq for NGC 4889, M, = 9.7l^;^ x 10^ M© 



20 



S 10' 



IS 

a 



10^ 



10* 



N4889 



N3842 

I S--?^ 



I — €& iSf 1 ^ 



-^^ 



• 



'^r^' 
• 
• 



^2832 ^ 
N7768 • 



/„\ McConnell 2011b •ii^ 

V**/ Belfiori2012 Q 



10" t N3842 H— ^^ 



10" 



10' 



N4889 H— (H?H 1^ 



• 



1^ 



1 L N^^-'^ 



N7768 H- 



tf- 



1 1 



(b) 



McConnell 2011b -Cf 
Gultekin 2009a Q 



■ 10'"- 



10" 



10' 



N4889 



I G II 1?-^-i 



h 



-^ 



N3842 



? 



"=r 



f-7^' — f 






N7768 ^ 



(c) 



Magorrian 1998 -if 
Beifiori2012 Q 



10" 



10" 



10' 



10" 



10" 



10' 



10" 



10" 



10' 



Predicted Mjjh from a (M^^ J 



Predicted Mbh from Lv (M^J Predicted Mbh from M,,^,, (M^J 
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for NGC 3842, M, ^ LStoJ ^ 10^ A/q for NGC 7768, 
and M, < 9.0 x lO^Af© for NGC 2832. The orbit mod- 
els also determine the stellar mass-to-light ratio Mi,/ Lb,, 
which we convert to Mi,/ Ly as described in Section [ 



We find Mi,/LR = 5.9 ± 1.7 MgL"^^ {MjLy = 7.4 ± 
2.1 MeL-^y) for NGC 4889, Mi,/LR = 5.2±0.8 MqL'^j^ 
{Mi,/Lv = 7.1 ± 1.1 M^L-^y) for NGC 3842, A/^/L 



^R 



5.7 ± l.lAfoL-;^^ {Mi,/Lv = 8.8 ± 1.7AfoL-V) for 
NGC 7768, and Mi,/ Lr = 7. 6^0^^ MqL'^j^ (Af*/Ly = 



9.7±1.0A^mL 



Q^Q,V. 



for NGC 2832. 



The orbital distributions from the best-fitting models 
indicate that NGC 4889 and NGC 3842 are depleted of 
radial orbits, out to radii of approximately 1 kpc. These 
extended tangential bias regions have similar sizes to 
the galaxies' photometric cores, and may be linked to 
the extremely massive black hole {M, ~ 10^" Mq) in 
each galaxy. NGC 7768 has a black hole mass of only 
~ 10^ Mq and a correspondingly modest tangential bias 
region. 

Several correlations between black hole mass and scalar 
galaxy properties have been explored by numerous au- 
thors. Below we focus on the Af, — a relation, the V- 
band M, — L relation, and the relation between Af, and 
bulge stellar mass, A/^. We estimate Mi, for each BCG 
by multiplying Ly by Mi,/ Ly- 

In Figure 1181 we compare our measured values of A/, 
in NGC 4889, NGC 3842, and NGC 7768, and our up- 
per limit for NGC 2832, to predicted black hole masses 
from the M, — a, M, — Ly, and M, — Mi, relations. 
In pa rticular, we display the a - and L-based predictions 
from iMcConnell et alj ()2011b[ ). who use the most up- 
to-date sample of directly measured black hole masses. 
The scaling relations derived from this sample are steep- 



ened by the large black hole masses in NGC 4889 and 
NGC 3842, as weU as upward revisions of M, in M87 
and M60 (Gebhardt & Thomas 2009; Shcn & Gcbhardl 
[20T0I: [Gebhardt et'aLllSoni) . The M, - a relation is 
additionally steepened by the inclusion of eight late- 
type galaxies with ma s er-based mea. s ureni e nts of M, 
(iKontradko et al.l [20081 'Greene et a l.l [2OIOI: iKuq et alJ 
120 lit) . As a result, the McConncU et al.l (|2011bD rela- 
tions predict the largest values of Af, in BCGs. For 
comparison, we have considered a large number of re- 
ported power-law fits to earlier galaxy samples and have 
determined which fits yield the lowest predicted values 
of M, for our BCGs. These minimum predictions are 
also displayed in Figure [181 in order to illustrate the full 
range of investigations to date. For the Af, — Mi, re- 
lation, we display the predi cted black hol e mass es from 
IMaeorrian et alJ ()1998f ) and IBeifiori et al.l (|2012D . which 
span the full range of Af*-based predictions. Along with 
the Af, values predicted from the mean relations, we con- 
sider the intrinsic scatter, eoj in log(Af,) at fixed a, Ly, 
or Mi,. The horizontal error bars in Figure [TSj essentially 
illustrate eg; measurement errors in a, Ly, and Af* have 
much smaller effects on the predicted black hole masses. 

The left panel of Figure \TE\ indicates that our measure- 
ments of Af, in NGC 4889, NGC 3842, and NGC 7768 are 
all greater than the predicted values from various fits to 
the Af, - CT relation. For NGC 4889 and NGC 3842, the 
discrepancy exceeds the intrinsic scatter in Af, regard- 
less of which fit to Af, - cr we select. For jMGC 7768, Af. 
exceeds the mean A^« — a rela tions from IGiiltekinet alJ 
(|2009al) . iGraham et all (|201l[ ) and IBeifiori et al.l (|2012f ) 
by more than the corresponding estimates of eoj but is 
consi stent with the pow er-law fit and eg estimate from 
iMcConnehet al.l (|2011bl) . 

The predicted black hole masses from M, — Ly, in- 
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dicated in the middle panel of Figure I18[ are uniformly 
more massive than the c-based predictions. Our large 
measurement errors for M, in NGC 4889, combined 
with estimates of cq, yield agreement between the mea- 
sured black hole mass and the M, — Ly relation of 
iMcConnell et alJ (|2011b[ ). However, the black hole in 
NGC 3842 is significantly more massive than predicted 
by any version of the M, — Ly relation. For NGC 7768, 
our measurement of M, agrees with all versions of the 
M, — Ly relation. 

The right panel of Figure [18] shows that the M, — M^, 
relation offers the widest range of predicted black hole 
mass es from the literature. Only the historical relation 
from lMagorrian et al.l (|1998[ ) is consistent with the mea- 
sured black hole masses of both NGC 4889 and NGC 
3842. Mor e recent estim ates of M, — Mi, by iSani et aD 
(|2011h and iBeifiori et al.l ()2012l ). however, predict black 
hole masses an order of magnitude lower and are consis- 
tent with our measurement in NGC 7768. 

Our upper limit for M, in NGC 2832 is consistent with 
the full range of predicted black hole masses from a, Ly, 
and Mi,. 

There has been recent dispute over the existence of 
a fundamental correlation between Af, and the Virial 
mass or circular velocity of the host galaxy's dark mat- 
ter halo (e.g., Ferrarese 2002; Kormendy & Bender 2011; 
Volonteri et al. 2011; also Beifiori et al. 2012 and refer- 
ences therein). We can compare NGC 4889 and NGC 
3842 to this correlation by using Vc from their best- 
fitting LOG dark matter halos, or by computing the 
Virial mass of the best-fitting NFW halo for NGC 3842 
(M200 = 8.5 X 10^^ Mq). Using either metric, we find 
that the black holes in NGC 4889 and NGC 3842 are 
several times more massive than those in galaxies with 
comparable dark matter halos. 

The present-day black hole mass function is a testable 
prediction of galaxy evolution and feedback models. In 
particular, different models make divergent predictions 
for the number de nsity of b lack holes with M, ~ 10^ — 
10^° M© (see, e.g.. |Kelly~fc M erloni 2012, and references 
therein). Our measurements of M, ~ 10^" M© in NGC 
4889 and NGC 3842, plus the 6.3 x lO^M© black hole 
in M87 (IGebhardt et all [Ml and th e 4.7 x IO^Mq 
black hole in M60 (jShen fc Gebhardtl I2010D . place a 
lower limit of four galaxies with M, > 10^'^ M© within 
a local volume ~ 10^ Mpc'^. This is slightly higher 
than the space density predicted from models of propor- 
tional bla ck hole and spheroid growth in major merg- 
ers (e.g., IHopkins et all I2008D . Several other BCGs 
lie within 100 Mpc, and the discovery of more black 
holes with M, > 10^'^ M© could a pproac h the pre- 
dicted space density from models by iShenI (|2009f ). in 
which all quasars follow a universal light curve and the 
most massive black ho l es acc rete gas until redshifts ~ 1. 
INataraian fc Treisteri (|2009l ) have predicted the black 
hole mass function from the X-ray luminosity function 
of active galactic nuclei (AGN). Our lower limit is con- 
sistent with their prediction using t he original luminos- 
ity function from lUeda et all (|2003( ) but does not sup- 
port their hypothesis of a truncated lu minosity function 
for AGN with Af. > 10^. Models bv lYoo et all (pOQTf ) 
find that black hole-black hole mergers can increase Af, 
by factors ^ 2 in galaxy clusters, with the largest black 



holes reaching 1 — 1.5 x 10-'^° Af©. Black hole mergers also 
dominate the growth of the larg est black holes in recent 
models bv lFanidakis et all (|2QllD . but their prescriptions 
for quasar-mode (near-Eddington) and radio-mode (low- 
Eddington) accretion predict relatively low space densi- 
ties for black holes with Af, > 10^ Af©. 

The steep L — a relation for BCGs may arise from 
extreme size evolu tion, either via low-angular mom en- 
tum mergers fe.g.. IBovlan-Kolchin et al]|2005l I2OO60 or 
via can nibalis m of numero u s low-den sity sys tems (e.g. , 
Ostriker fc Tremainel 119751: lOstriker fc HausmanI 119771: 



Hopkins et al.ll2010l: lOser et al.ll2010[ ). In both scenar 
ios, size evolution is driven by gas-poor or "dry" pro 
cesses. Dry major mergers should preserve any initial 
relation between Af, and stellar mass. For black hole 
growth to mirror stellar mass growth in the cannibal- 
ism scenario, the devoured satellites must indeed contain 
central black holes, and dynamical friction must deliver 
them efficiently to the center of the BCG. The forms of 
the Af, — a, Af, — L, and Af, — Af^, relations for BCGs 
could provide additional insight to their progenitors and 
growth mechanisms. 

Our BCGs are not the only systems that have 
been observed to deviate from the black hole scal- 
ing relations. The Af, — a, Af, — L, and Af, — Af^, 
relations each have outliers that span a variety of 
galaxy masses and environments (e.g., Nowak et al. 2008 
Greene et al. 2010; Nowak et al. 2010; Rcincs ct al. 2011 
iBogdan et al.|[20T2[ ). In fact, some merging models pre 
diet that the most massive galaxies will exhibit smaller 
deviations from the mean black hole scaling relations, as 
a consequence of experiencing gre ater numbers of hierar- 
chical mergers ( Peng 2007; Volo nteri &: NataraianI [20091 : 
iJahnke fc Macciol 120111) . In this scenario, black hole 
growth and galaxy growth need not be coupled through 
direct feedback processes. On the other hand, systematic 
offsets for Af, in BCGs could point to unique evolution- 
ary processes near cluster centers. A larger sample of 
BCGs with measured Af, is necessary to distinguish be- 
tween systematic trends and random scatter. 

We aim to use high-resolution and wide-field kinematic 
data of several more galaxies to explore the high-mass 
slopes of the empirical correlations between galaxies and 
supermassive black holes, and to quantify the total and 
intrinsic scatter in Af, for massive galaxies in different 
environments. This investigation and future studies are 
important tests for the hypothesis of universal scaling 
relations between black holes and their host galaxies, and 
for galaxy evolution models using different modes and 
timescales for black hole accretion. 
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